Molecular insights into the variability and pathogenicity of Fusarium odoratissimum, the causal agent of Panama wilt disease in banana.

Fusarium wilt or Panama disease of banana caused by the hemibiotroph fungus, Fusarium odoratissimum, also known as F. oxysporum f.sp. cubense Tropical Race 4 is a serious threat to banana production worldwide. Being the world's largest grower and the origins of bananas in its northeast region, India is particularly vulnerable to this deadly fungus. In the present study, a total of 163 Fusarium isolates from infected banana were characterized for their pathogenic traits. Considering the variability in the Fusarium, the contaminated banana plants were collected from five districts of Uttar Pradesh and Bihar, two major primary infection states of India. All the isolates were screened using universal and specific primers to identify the F. odoratissimum strains. The identified F. odoratissimum strains were subjected to in vivo pathogenicity assessment using the susceptible banana cultivar 'Grand Naine'. The identified six most virulent strains were further characterized for their pathogenicity via in vivo bipartite interaction in terms of biochemical assays. Assessment of in vivo pathogenicity through qRT-PCR for three pathogenesis responsive genes, Six 1a (Secreted in xylem), Snf (Sucrose non-fermenting) and ChsV (Chitinase V), ascertained that the identified F. odoratissimum strains exhibit both intra- and inter-specific variability. The variability of F. odoratissimum strains signifies its importance for the assessment of spread of infection at specific sites to enable efficient management strategy of Fusarium wilt in banana.

Genome-wide association study of fiber yield-related traits uncovers the novel genomic regions and candidate genes in Indian upland cotton (Gossypium hirsutum L.).

Upland cotton (Gossypium hirsutum L.) is a major fiber crop that is cultivated worldwide and has significant economic importance. India harbors the largest area for cotton cultivation, but its fiber yield is still compromised and ranks 22nd in terms of productivity. Genetic improvement of cotton fiber yield traits is one of the major goals of cotton breeding, but the understanding of the genetic architecture underlying cotton fiber yield traits remains limited and unclear. To better decipher the genetic variation associated with fiber yield traits, we conducted a comprehensive genome-wide association mapping study using 117 Indian cotton germplasm for six yield-related traits. To accomplish this, we generated 2,41,086 high-quality single nucleotide polymorphism (SNP) markers using genotyping-by-sequencing (GBS) methods. Population structure, PCA, kinship, and phylogenetic analyses divided the germplasm into two sub-populations, showing weak relatedness among the germplasms. Through association analysis, 205 SNPs and 134 QTLs were identified to be significantly associated with the six fiber yield traits. In total, 39 novel QTLs were identified in the current study, whereas 95 QTLs overlapped with existing public domain data in a comparative analysis. Eight QTLs, qGhBN_SCY_D6-1, qGhBN_SCY_D6-2, qGhBN_SCY_D6-3, qGhSI_LI_A5, qGhLI_SI_A13, qGhLI_SI_D9, qGhBW_SCY_A10, and qGhLP_BN_A8 were identified. Gene annotation of these fiber yield QTLs revealed 2,509 unique genes. These genes were predominantly enriched for different biological processes, such as plant cell wall synthesis, nutrient metabolism, and vegetative growth development in the gene ontology (GO) enrichment study. Furthermore, gene expression analysis using RNAseq data from 12 diverse cotton tissues identified 40 candidate genes (23 stable and 17 novel genes) to be transcriptionally active in different stages of fiber, ovule, and seed development. These findings have revealed a rich tapestry of genetic elements, including SNPs, QTLs, and candidate genes, and may have a high potential for improving fiber yield in future breeding programs for Indian cotton.

De Novo Hybrid Assembled Draft Genome of Commiphora wightii (Arnott) Bhandari Reveals Key Enzymes Involved in Phytosterol Biosynthesis.

Genome sequence and identification of specific genes involved in the targeted secondary metabolite biosynthesis are two essential requirements for the improvement of any medicinal plant. Commiphora wightii (Arnott) Bhandari (family: Burseraceae), a medicinal plant native to Western India, produces a phytosterol guggulsterone, which is useful for treating atherosclerosis, arthritis, high cholesterol, acne, and obesity. For enhanced guggulsterone yield, key genes involved in its biosynthesis pathway need to be predicted, for which the genome sequence of the species is a pre-requisite. Therefore, we assembled the first-ever hybrid draft genome of C. wightii with a genome size of 1.03 Gb and 107,221 contigs using Illumina and PacBio platforms. The N50 and L50 values in this assembled genome were ~74 Kb and 3486 bp, respectively with a guanine-cytosine (GC) content of 35.6% and 98.7%. The Benchmarking Universal Single Copy Ortholog (BUSCO) value indicated good integrity of assembly. Analysis predicted the presence of 31,187 genes and 342.35 Mb repeat elements in the genome. The comparative genome analysis of C. wightii with relevant orthogroups predicted a few key genes associated with phytosterol biosynthesis and secondary metabolism pathways. The assembled draft genome and the predicted genes should help the future variety development program with improved guggulsterone contents in C. wightii.

Chloroplast Genome of Lithocarpus dealbatus (Hook.f. & Thomson ex Miq.) Rehder Establishes Monophyletic Origin of the Species and Reveals Mutational Hotspots with Taxon Delimitation Potential

There is phylogenetic ambiguity in the genus Lithocarpus and subfamily Quercoideae (Family: Fagaceae). Lithocarpus dealbatus, an ecologically important tree, is the dominant species among the Quercoideae in India. Although several studies have been conducted on the species' regeneration and ecological and economic significance, limited information is available on its phylo-genomics. To resolve the phylogeny in Quercoideae, we sequenced and assembled the 161,476 bp chloroplast genome of L. dealbatus, which has a large single-copy section of 90,732 bp and a small single-copy region of 18,987 bp, separated by a pair of inverted repeat regions of 25,879 bp. The chloroplast genome contained 133 genes, of which 86 were protein-coding genes, 39 were transfer RNAs, and eight were ribosomal RNAs. Analysis of repeat elements and RNA editing sites revealed interspecific similarities within the Lithocarpus genus. DNA diversity analysis identified five highly diverged coding and noncoding hotspot regions in the four genera, which can be used as polymorphic markers for species/taxon delimitation across the four genera of Quercoideae viz., Lithocarpus, Quercus, Castanea, and Castanopsis. The chloroplast-based phylogenetic analysis among the Quercoideae established a monophyletic origin of Lithocarpus, and a closer evolutionary lineage with a few Quercus species. Besides providing insights into the chloroplast genome architecture of L. dealbatus, the study identified five mutational hotspots having high taxon-delimitation potential across four genera of Quercoideae.

Identification of Stable and Multiple Environment Interaction QTLs and Candidate Genes for Fiber Productive Traits Under Irrigated and Water Stress Conditions Using Intraspecific RILs of Gossypium hirsutum var. MCU5 X TCH1218.

Cotton productivity under water-stressed conditions is controlled by multiple quantitative trait loci (QTL). Enhancement of these productivity traits under water deficit stress is crucial for the genetic improvement of upland cotton, Gossypium hirsutum. In the present study, we constructed a genetic map with 504 single nucleotide polymorphisms (SNPs) covering a total span length of 4,416 cM with an average inter-marker distance of 8.76 cM. A total of 181 intra-specific recombinant inbred lines (RILs) were derived from a cross between G. hirsutum var. MCU5 and TCH1218 were used. Although 2,457 polymorphic SNPs were detected between the parents using the CottonSNP50K assay, only 504 SNPs were found to be useful for the construction of the genetic map. In the SNP genotyping, a large number of SNPs showed either >20% missing data, duplication, or segregation distortion. However, the mapped SNPs of this study showed collinearity with the physical map of the reference genome (G. hirsutum var.TM-1), indicating that there was no chromosomal rearrangement within the studied mapping population. RILs were evaluated under multi-environments and seasons for which the phenotypic data were acquired. A total of 53 QTL controlling plant height (PH), number of sympodial branches, boll number (BN), and boll weight (BW) were dissected by QTL analysis under irrigated and water stress conditions. Additionally, it was found that nine QTL hot spots not only co-localized for more than one investigated trait but were also stable with major QTL, i.e., with > 10% of phenotypic variation. One QTL hotspot on chromosome 22 flanked by AX-182254626-AX-182264770 with a span length of 89.4 cM co-localized with seven major and stable QTL linked to a number of sympodial branches both under irrigated and water stress conditions. In addition, putative candidate genes associated with water stress in the QTL hotspots were identified. Besides, few QTL from the hotspots were previously reported across various genetic architects in cotton validating the potential applications of these identified QTL for cotton breeding and improvement. Thus, the major and stable QTL identified in the present study would improve the cotton productivity under water-limited environments through marker-assisted selection.

GBS-SNP and SSR based genetic mapping and QTL analysis for drought tolerance in upland cotton.

A recombinant inbred line mapping population of intra-species upland cotton was generated from a cross between the drought-tolerant female parent (AS2) and the susceptible male parent (MCU13). A linkage map was constructed deploying 1,116 GBS-based SNPs and public domain-based 782 SSRs spanning a total genetic distance of 28,083.03 cM with an average chromosomal span length of 1,080.12 cM with inter-marker distance of 10.19 cM.A total of 19 quantitative trait loci (QTLs) were identified in nine chromosomes for field drought tolerance traits. Chromosomes 3 and 8 harbored important drought tolerant QTLs for chlorophyll stability index trait while for relative water content trait, three QTLs on chromosome 8 and one QTL each on chromosome 4, 12 were identified. One QTL on each chromosome 8, 5, and 7, and two QTLs on chromosome 15 linking to proline content were identified. For the nitrate reductase activity trait, two QTLs were identified on chromosome 3 and one on each chromosome 8, 13, and 26. To complement our QTL study, a meta-analysis was conducted along with the public domain database and resulted in a consensus map for chromosome 8. Under field drought stress, chromosome 8 harbored a drought tolerance QTL hotspot with two in-house QTLs for chlorophyll stability index (qCSI01, qCSI02) and three public domain QTLs (qLP.FDT_1, qLP.FDT_2, qCC.ST_3). Identified QTL hotspot on chromosome 8 could play a crucial role in exploring abiotic stress-associated genes/alleles for drought trait improvement. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-021-01041-y.

Genetic diversity and population structure of Arabidopsis thaliana along an altitudinal gradient.

The natural genetic variation within a plant species is primarily a consequence of its phylogeography and evolutionary history. This variation largely determines its present-day population structure. Arabidopsis thaliana, as a model plant, has been studied in great detail including its probable origin, local as well as global genetic diversity pattern, population structure, adaptation, etc. However, no such studies have so far been reported from the Indian Himalayan region. Here, we describe a comprehensive study on the genetic diversity and population structure of A. thaliana from an altitudinal range of 700-3400 m above mean sea level the highest altitudinal range reported so far. We also compare these populations with previously reported worldwide populations. A total of 48 accessions representing six populations were analysed using 19 microsatellites and 11 chloroplast markers. Genetic diversity analysis indicated populations to be highly diverse and comparable with worldwide populations. STRUCTURE, principal coordinate and isolation by distance (IBD) analyses showed that genetic variation in different populations is structured at geographical and altitudinal level. Further analyses indicate that these populations are genetically distinct from the rest of the world populations. Different parameters of the demographic expansion model support a rapid expansion. Based on mismatch distribution, the initial time of expansion of west Himalayan populations was found to be about 130 000 years. Bayesian analysis of divergence time indicated that these populations have a long evolutionary history in this region. Based on the results of genetic diversity parameters, demographic expansion and divergence time estimation, it appears that west Himalayan populations may be the source of the west-east expansion model.

Microarray-based large scale detection of single feature polymorphism in Gossypium hirsutum L.

Microarrays offer an opportunity to explore the functional sequence polymorphism among different cultivars of many crop plants. The Affymetrix microarray expression data of five genotypes of Gossypium hirsutum L. at six different fibre developmental stages was used to identify single feature polymorphisms (SFPs). The background corrected and quantile-normalized log2 intensity values of all probes of triplicate data of each cotton variety were subjected to SFPs call by using SAM procedure in R language software. We detected a total of 37,473 SFPs among six pair genotype combinations of two superior (JKC777 and JKC725) and three inferior (JKC703, JKC737 and JKC783) using the expression data. The 224 SFPs covering 51 genes were randomly selected from the dataset of all six fibre developmental stages of JKC777 and JKC703 for validation by sequencing on a capillary sequencer. Of these 224 SFPs, 132 were found to be polymorphic and 92 monomorphic which indicate that the SFP prediction from the expression data in the present study confirmed a ~58.92% of true SFPs. We further identified that most of the SFPs are associated with genes involved in fatty acid, flavonoid, auxin biosynthesis etc. indicating that these pathways significantly involved in fibre development.

Large-scale resource development in Gossypium hirsutum L. by 454 sequencing of genic-enriched libraries from six diverse genotypes.

The sequence information has been proved to be an essential genomic resource in case of crop plants for their genetic improvement and better utilization by humans. To dissect the Gossypium hirsutum genome for large-scale development of genomic resources, we adopted hypomethylated restriction-based genomic enrichment strategy to sequence six diverse genotypes. Approximately 5.2-Gb data (more than 18.36 million reads) was generated which, after assembly, represents nearly 1.27-Gb genomic sequences. We predicted a total of 93,363 gene models (21,399 full length) and identified 35,923 gene models which were validated against already sequenced plant genomes. A total of 1,093 transcription factor-encoding genes, 3,135 promoter sequences and 78 miRNA (including 17 newly identified in Gossypium) were predicted. We identified significant no. of molecular markers including 47,093 novel simple sequence repeats and 66,364 novel single nucleotide polymorphisms. In addition, we developed NBRI-Comprehensive Cotton Genomics database, a web resource to provide access of cotton-related genomic resources developed at NBRI. This study contributes considerable amount of genomic resources and suggests a potential role of genic-enriched sequencing in genomic resource development for orphan crop plants.

Development and characterization of genomic and expressed SSRs for levant cotton (Gossypium herbaceum L.).

Four microsatellite-enriched genomic libraries for CA(15), GA(15), AAG(8) and ATG(8) repeats and transcriptome sequences of five cDNA libraries of Gossypium herbaceum were explored to develop simple sequence repeat (SSR) markers. A total of 428 unique clones from repeat enriched genomic libraries were mined for 584 genomic SSRs (gSSRs). In addition, 99,780 unigenes from transcriptome sequencing were explored for 8,900 SSR containing sequences with 12,471 expressed SSRs. The present study adds 1,970 expressed SSRs and 263 gSSRs to the public domain for the use of genetic studies of cotton. When 150 gSSRs and 50 expressed SSRs were tested on a panel of four species of cotton, 68 gSSRs and 12 expressed SSRs revealed polymorphism. These 200 SSRs were further deployed on 15 genotypes of levant cotton for the genetic diversity assessment. This is the first report on the successful use of repeat enriched genomic library and expressed sequence database for microsatellite markers development in G. herbaceum.